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Abstract. Stars form by gravoturbulent fragmentation of interstellar gas clouds. The supersonic 
turbulence ubiquitously observed in Galactic molecular gas generates strong density fluctuations 
with gravity taking over in the densest and most massive regions. Collapse sets in to build up 
stars and star clusters. 

Turbulence plays a dual role. On global scales it provides support, while at the same time 
it can promote local collapse. Stellar birth is thus intimately linked to the dynamic behavior 
of parental gas clouds, which governs when and where protostellar cores form, and how they 
contract and grow in mass via accretion from the surrounding cloud material to build up stars. 
The equation of state plays a pivotal role in the fragmentation process. Under typical cloud 
conditions, massive stars form as part of dense clusters following the "normal" mass function 
observed, e.g. in the solar neighborhood. However, for gas with an effective polytropic index 
greater than unity star formation becomes biased towards isolated massive stars. This is relevant 
for understanding the properties of zero-metallicity stars (Population III) or stars that form 
under extreme environmental conditions like in the Galactic center or in luminous starbursts. 
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1. Introduction 

Identifying the physical processes that determine the masses of stars and their sta- 
tistical distribution, usually known as the initial mass function (IMF), is a fundamental 
problem in star formation research. Observations in the solar vicinity and in nearby young 
star clusters suggest the existence of a characteristic stellar mass M cri . The IMF peaks at 
this characteristic mass which is typically a few tenths of a solar mass. For larger masses 
the IMF has a nearly power-law form and it declines rapidly towards smaller masses 
(Scalo 1998; Kroupa 2002; Chabrier 2003). 

Although the IMF in our Milky Way has been derived from vastly different regions, 
ranging from the immediate neighborhood of our Sun to dense and distant stellar clusters, 
its key features seem to be strikingly universal to all determinations (e.g. Kroupa 2001, 
Meyer et al. 2004 and references therein). However, the environmental conditions in the 
considered regions vary enormously. If the formation of stars and the resulting mass 
spectrum strongly depend on initial conditions, there would thus be no reason for the 
IMF to be universal. A derivation of the characteristic stellar mass that is based on 
fundamental atomic and molecular physics would therefore be highly desirable. 

On the other hand, there are hints for IMF variations in more extreme environments 
and/or at different metallicities. For example, HST/WFPC2 observations of the LMC 
by Gouliermis, Brander & Henning (2004) indicate that local conditions seem to favor 
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the formation of higher mass stars (top-heavy IMF) in associations, and not in the 
background field. Similar conclusions may be drawn for the stellar mass spectrum in 
the Galactic Center (Mark Morris, this proceedings; or Stolte et al. 2003 for the Arches 
Cluster) or in starburst galaxies (e.g. Elmegreen 2004, and references therein). 

There are various ways to approach the formation of stars and star clusters from a 
theoretical point of view. Currently, the most successful models are based on the in- 
terplay between self-gravity and the interstellar turbulence ubiquiteously observed in 
Galactic molecular clouds. In this picture, stars form by a process we call gravoturbulent 
fragmentation (see, e.g. the reviews by Larson 2003, or Mac Low & Klessen 2004). How- 
ever, despite a wealth of observational data for the Milky Way (e.g. Fuller & Myers 1992; 
Goodman et al. 1998), the nature of interstellar turbulence and the thermodynamic state 
of the star-forming gas (as characterized by the equation of state - EOS) remain major 
theoretical problems in understanding the stability and collapse of molecular clouds. The 
stiffness of the EOS can be largely responsible for the resulting density probability func- 
tion of interstellar gas in the turbulent ISM. In particular, the value of the polytropic 
index 7 when adopting an EOS of the form P oc p 7 strongly influences the compress- 
ibility of density condensations as well as the amount of clump fragmentation. Hence, 
the EOS is directly related to the IMF (Vazquez-Semadeni et al. 1996; Li, Klessen, & 
Mac Low 2003; Jappsen et al. 2005). In Spaans & Silk (2000, 2005) the properties of 
a polytropic EOS were investigated and it was found that the stiffness of the EOS de- 
pends strongly on the ambient metallicity and the infrared background radiation field 
produced by irradiated dust grains. It thus may vary considerably in different galactic 
environments. 

The structure of this paper is as follows: In the next section, SJ2 we briefly review the 
basic concepts of gravoturbulent fragmentation and discuss its relation to stellar birth. In 
JHwe consider the fragmentation behavior of self-gravitating, supersonically turbulent gas 
with a polytropic EOS. Small values of the polytropic exponent 7 lead to the formation 
of dense clusters of low-mass stars, while 7 > 1 introduces a bias towards more massive 
and isolated stars. In 21 we discuss the thermal properties of the star-forming gas in our 
Galaxy and their implications for the IMF in the solar neighborhood. Finally, in we 
present first numerical calculations of gravoturbulent fragmentation in starburst galaxies 
suggesting that the IMF in very extreme environments may be top-heavy. 

2. The Concept of Gravoturbulent Fragmentation 

Supersonic turbulence plays a dual role in star formation. While it usually is strong 
enough to counterbalance gravity on global scales it will usually provoke collapse locally. 
This process is discussed in detail in the reviews by Larson (2003) and Mac Low & 
Klessen (2004). Turbulence establishes a complex network of interacting shocks, where 
regions of high-density build up at the stagnation points of convergent flows. These gas 
clumps can be dense and massive enough to become gravitationally unstable and collapse 
when the local Jeans length becomes smaller than the size of the fluctuation. However, 
the fluctuations in turbulent velocity fields are highly transient. They can disperse again 
once the converging flow fades away (Vazquez-Semadeni et al. 2005). Even clumps that 
are strongly dominated by gravity may get disrupted by the passage of a new shock front 
(Mac Low et al. 1994). 

For local collapse to result in the formation of stars, Jeans unstable, shock-generated, 
density fluctuations must collapse to sufficiently high densities on time scales shorter 
than the typical time interval between two successive shock passages. Only then do they 
'decouple' from the ambient flow pattern and survive subsequent shock interactions. The 
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shorter the time between shock passages, the less likely these fluctuations are to survive. 
The overall efficiency of star formation depends strongly on the wavelength and strength 
of the driving source (Klessen et al. 2000; Heitsch et al. 2001). Both regulate the amount 
of gas available for collapse on the sonic scale where turbulence turns from supersonic to 
subsonic (Vazquez-Semadeni et al. 2003). 

The velocity field of long-wavelength turbulence is dominated by large-scale shocks 
which are very efficient in sweeping up molecular cloud material, thus creating massive 
coherent structures. These exceed the critical mass for gravitational collapse by far. The 
situation is similar to localized turbulent decay, and quickly a cluster of protostellar 
cores builds up (e.g. Klessen et al. 1998; Klessen 2001; Bate et al. 2003; Li et al. 2004; or 
Bonnell, this proceedings). A prominent example is the Trapezium Cluster in Orion with 
a few thousand young stars (Hillenbrand 1997). However, this scenario also applies to the 
Taurus star-forming region (Hartmann 2002) which is historically considered as a case of 
isolated stellar birth. Its stars have formed almost simultaneously within several coherent 
filaments which apparently are created by external compression (see Ballesteros-Paredes 
et al. 1999). This renders it a clustered star- forming region in the sense of the above 
definition. 

The efficiency of turbulent fragmentation is reduced if the driving wavelength decreases. 
There is less mass at the sonic scale and the network of interacting shocks is very tightly 
knit. Protostellar cores form independently of each other at random locations throughout 
the cloud and at random times. There are no coherent structures with multiple Jeans 
masses. Individual shock generated clumps are of low mass and the time interval between 
two shock passages through the same point in space is small. Hence, collapsing cores are 
easily destroyed again. Altogether star formation is inefficient, and stars are dispersed 
throughout the cloud. 

3. Gravoturbulent Fragmentation in Polytropic Gas 

As indicated above, the thermodynamic properties of the star-forming gas may play a 
crucial role in determining the stellar mass spectrum. As a first and crude approximation, 
the balance between heating and cooling in a molecular cloud can be described by a 
polytropic EOS, P = Kp 1 ', where if is a constant, and P, p and 7 are thermal pressure, 
gas density, and polytropic exponent, respectively. Depending on the environment, the 
exponent 7 can cover the range 0.2 < 7 < 1.4 in the interstellar medium (e.g. Spaans & 
Si lk 2000 2 005; S calo & Biswas 2002). 

iLi et al 1 ll2003h carried out detailed numerical calculations to determine the effects of 
different EOS's on gravoturbulent fragmentation by varying 7 in steps of 0.1 in otherwise 
identical simulations. Figure ^ illustrates how a low exponent 7 leads to the build-up of 
a dense cluster of stars, while high values of 7 result in isolated star formation. It also 
shows that the spectra of both the gas clumps and protostars change with 7. In low-7 
models, the mass distribution of the collapsed protostellar cores at the high-mass end is 
roughly log-normal. As 7 increases, fewer but more massive cores emerge. When 7 > 1.0, 
the distribution is dominated by high mass protostars only, and the spectrum tends to 
flatten out. It is no longer described by either a log-normal or a power-law. The clump 
mass spectra, on the other hand, do show power-law behavior at the high mass side, even 
for 7 > 1.0. 

This suggests that in a low-7 environment stars tend to form in clusters and with small 
masses. On the other hand, massive stars can form in small groups or in isolation in gas 
with 7 > 1.0. 

The usefulness of theoretical predictions for the formation of isolated massive stars may 
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Figure 1. Top: 3-D distribution of the gas and protostars for different 7. Bottom: Mass spectra 
of gas clumps (thin lines) and of protostars (collapsed cores: hatched thick-lined histograms) for 
the corresponding cube above. The percentage shows the fraction of total mass accreted onto 
protostars. The vertical line shows the numerical resolution limit. Shown also are two p ower-law 
spectr a with v = —1.5 (dashed-line) and v — —2.33 (dotted line). Figure adopted from lLi et alJ 
(2003). 



be debated, since u sually, massive stars are found in clusters. It is therefore interesting to 
note, that recently [Lamers et al.l (|2002t) reported observations of i solated massive stars 
or very small groups of massive stars in the bulge of M51. Also iMassevI |.2002) finds 
massive, apparently isolated field stars in both the Large and Small Magellanic Clouds. 
This is consistent with our models assuming 7 > 1.0. 



4. Thermal Properties of Galactic Gas and Implications for the IMF 

Observational evidence indicates tha t dense prestellar cloud cores show a rough balance 
between gravity and thermal pressure l|Benson fc Mver 

3 ll 98fiUMvers et ;l lJll«)?)lh. Thus, 
the thermodynamic properties of the gas play an important role in determining how 
dense star- forming regions in molecular clouds collapse and fragment. Observational and 
theoretical studies of the thermal properties of collapsing clouds both indicate that at 
densities below 10 -18 gcm~ 3 , roughly corresponding to a number density of n — 2.5 x 
10 5 cm -3 , the temperature decreases with increa sing density. This is due to the strong 
dependence of molecular cooling rates on density ijKovama k. Inutsuk a 2000). Therefore, 
the polytropic exponent 7 is below unity in this density regime. At densities above 
10~ 18 gcm~ 3 , the gas becomes thermally coupled to the dust grains, which then control 
the temperature by far-infrared thermal emission. The balance between compressional 
heating and thermal cooling by dust causes the temperature to increase again slowly with 
increasing density. Thus th e temperature-dens i ty relation can be approxi mated with 7 
above unity in this regime l|LarsorJll985l l2005t ISpaans fc SilldkoOOt 12005^1 . Changing 7 
from a value below unity to a va lue above unity results in a minimum temperature at 
the critical density. As shown bv iLi et "ahl l|2003j) . gas fragments efficiently for 7 < 1.0 
and less efficiently for higher 7. Thus, the Jeans mass at the critical density defines a 
characteristic mass for fragmentation, which may be r elated to the peak of the IMF. 
This hypothesis has been tested bv I Jappsen et alJ J2005), who studied the effects of 
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Figure 2. Mass spectra of protostellar cores for four models with critical densities in the range 
4.3 x 10 4 cm" 3 < n c ^ 4.3 x 10 7 cm . We show two phases of evolution, when about 10% 
and 30% of the total mass has been accreted onto protostars. The vertical solid line s hows the 
posit ion of the median mass. The dotted line serves as a reference to the Salpeter value (ISalpeterl 
1955). The dashed line indicates the mass resolution limit. (Adopted from Jappsen et al. 2005.) 



using a piecewise polytropic EOS with a polytropic exponent that changes from 7 = 0.7 
to 7 = 1.1 at critical densities in the range 4.3 x 10 4 cm~ 3 ^ n c ^ 4.3 x 10 7 cm~ 3 . Figure 
illustrates how the characteristic mass M c h for fragmentation dep ends on the adopted 
critical density. Close st correspondence with the observed IMF (see, IScalolll998t iKroupal 
2002] IChabrier1l2003|) occurs for a critical density of 4.3 x 10 6 cm~ 3 and for stages of 
accretion around 30% and above. At high masses, the distribution of collapsed objects 
follows a Salpe ter-like power law. For comparison we indicate the Salpeter slope x ~ 1.3 
llSalpeteri ri955) where the IMF is defined by d N/ d log mocm x . For masses around the 
median mass the distribution exhibits a small plateau and then falls off towards smaller 
masses. The top-row model where the change in 7 occurs below the initial mean density, 
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shows a flat distribution with only a few, but massive cores. They reach masses up to 
10 M© and the minimal mass is about 0.3 M Q . The distribution becomes more peaked 
for higher n c and there is a shift to lower masses. 

Altogether, the characteristic mass M c h decreases with increasing critical density n c 
following the relation M c h oc n ( T - 5±01 . The density at which 7 changes from below unity 
to above unity thus selects a characteristic mass scale for fragmentation. For physical 
parameters appropriate for the local interstellar medium, the resulting mass spectrum 
exhibits striking similarities to the observed IMF. 

5. The IMF in Starburst Galaxies 

In this section, we present a preliminary and first approach to understanding the 
IMF in starburst galaxies using direct numerical simulations of gravoturbulcnt cloud 
fragmentation. 

Again, we employ a polytropic EOS, P oc p 1 1 where 7 is the polytropic index and p the 
mass density. We adopt a perfect gas equation of state, P oc pT g for the gas temperature 
T g , to write 7 as 

dlogT, 
cflogp 

This last step is justified (Scalo & Biswas 2002; Vazquez- Semadeni, Passot & Pouquet, 
1996) as long as the heating and cooling terms in the fluid energy equation can adjust to 
balance each other on a time-scale shorter than the time-scale of the gas dynamics (i.e., 
thermal equilibrium). The model described in Spaans & Silk (2000, 2005) is used in this 
work and the interested reader is referred to these papers for a detailed description of 
the various heating and cooling terms that influence the polytropic index. 

The EOS was computed for self-gravitating spherical clouds approximated as singular 
isothermal spheres (Neufeld, Lepp & Melnick 1995), for which the total hydrogen number 
density nn scales with column density Ah per unit of velocity as 

A H = 7.2 x 10 19 n° i 5 cm- 2 /(kms~ 1 ). 

For the starburst environments of interest here, we assume a cosmic ray ionization rate 
of 3 x 10~ 15 s _1 . Solar relative abundances are assumed (Asplund et al., 2004; Jenkins 
2004) while the overall metallicity is two times solar (Barthel 2004). We further adopt 
a MRN grain size distribution (Mathis, Rumpl & Nordsieck 1977) and assume that the 
dust abundance scales with the carbon abundance. For the starburst model AVtur = 3 
km/s is adopted to take the larger input of kinetic energy (e.g. through supernovae) into 
account. The dust temperature inside the model clouds is set by a fiducial background star 
formation rate of 100 M Q yr _1 which causes dust to be at temperatures of about T d = 70 
to Td = 30 K, depending on the amount of shielding. These values are consistent with 
dust temperatures determined by continuum observations of luminous infrared galaxies 
(Spinoglio et al., 2002; Klaas et al. 2001, 1997). No freeze-out of molecules is assumed. 
The latter assumption is justified because our dust is generally warmer than 20 K. 

An overabundance, relative to solar, of oxygen leads to a significantly stiffer EOS be- 
cause the opacity in oxygen bearing molecules like CO and water is higher. This increases 
the amount of line trapping in these opaque cooling lines. Furthermore, when the tem- 
perature of the dust is high (larger than 30 K), H 2 collisional de-excitation heating 
dominates the thermal budget because the water levels are pumped by far-infrared dust 
emission. This causes the gas heating to attain a n H density dependence (Takahashi, 
Hollenbach & Silk 1983; Spaans & Silk 2000, 2005). Consequently, one has a strong peak, 
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Figure 3. Starburst EOS adopted from Spaans & Silk (2005). 

7 w 1.4, around a density of ~ 10 4 cm" 3 . Since the dust grains are warm, gas-dust 
heating around 10 5 cm -3 causes the EOS to remain above unity, 7 w 1.1, beyond the 
H2O opacity peak. Note here that gas-dust heating/cooling scales as (T g — 7d)T° -5 njj. 
Hence, a starburst has a stiff EOS because of its high absolute metallicity and warm, 
Td ~ 30 — 70 K, dust. Figure [3] shows the resulting EOS and the fit that has been adopted 
for the hydrodynamic simulation. 

We adopt this EOS to model gravoturbulent star formation in a circumnuclear star- 
burst environment using smoothed particle hydrodynamics (SPH). This is a Lagrangian 
method to solve the equations of hydrodynamics, where the fluid is represented by an 
ensemble of particles, and flow quantities are obtained by averaging over an appropriate 
subset of SP H particles (Benz 199 0; Monaghan 1992). We use the parallel code GADGET, 
designed bv lSpringel et all ((2001) . The method is able to resolve large density contrasts 
as particles are free to move, and so the particle concentration increases naturally in high- 
density regio ns. To determine the re solution limit of our calculation we deploy a criterium 
proposed bv lBate fc BurkerU l)l997h . It requires the local Jeans mass at all densities to 
be represented by at least 100 particles, and is adequate for the problem considered here, 
where we follow the evolution of highly nonlinear density fluctuations created by super- 
sonic turbulen ce. Our version of the code fu rthermore replaces high-density cores with 
sink particles ijBate. Bonnell. & Pried Il995h . Sink particles can accrete gas from their 
surroundings while keeping track of mass and momentum. This enables us to follow the 
dynamic evolution of the system over many local free-fall timescales. 

We focus on a cubic molecular cloud region of 11.2 pc in size. The region contains 
80 000Mq of gas and has a mean particle density n — 1 000 cm" 3 . We drive turbulence 
uniformly on large scales, with wave numbers k in the range 1 ^ k 2 as described 
in Mac Low (1999). The energy input rate is adjusted to yield a constant turbulent 
Mach number A4 rms « 5. The particle number in the simulation presented here is N — 

8 000 000, which is thus one of the highest resolution star-formation calculation ever done 
with SPH. The critical density for sink particle formation is n c = 10 7 cm -3 . The mass 
of individual SPH particles is m = 0.01 M Q which is sufficient to resolve the minimum 
Jeans mass in the system Mj ss 1.5 Mq. Except for the EOS and th e particle number, 
the numerical set-up is identical to the study bv lJappsen et all fcOOfil) . 

As a preliminary result, we find that in the considered starburst galaxy, the mass 
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Figure 4. Mass spectrum of protostellar objects in a simulation focusing on a starburst envi- 
ronment. The size of the considered molecular cloud region is 11.2 pc and it contains 80 000 Mq. 
The system is depicted at early times when about 0.02% of the available gas is converted into 
stars. Comparison with Fig. [5] shows that the average mass (indicated by the vertical line) is 
shifted towards high values. The nominal resolution limit lies at IMq. 

spectrum of protostellar objects is biased towards high masses. This is illustrated in 
Fig. 01 The calculation thus supports the hypothesis that for extreme environmental 
conditions such as inferred for very dusty and IR-luminous circumnuclear starbursts the 
stellar initial mass function may indeed be top-heavy. 
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